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Phase separation and vortex states in binary mixture of Bose-Einstein condensates in 

the trapping potentials with displaced centers 
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The system of two simultaneously trapped codensates consisting of Rb atoms in two different 
hyperfine states is investigated theoretically in the case when the minima of the trapping potentials 
are displaced with respect to each other. It is shown that the small shift of the minima of the 
trapping potentials leads to the considerable displacement of the centers of mass of the condensates, 
in agreement with the experiment. It is also shown that the critical angular velocities of the vortex 
states of the system drastically depend on the shift and the relative number of particles in the 
condensates, and there is a possibility to exchange the vortex states between condensates by shifting 
the centers of the trapping potentials. 

PACS numbers: 03.75.Fi,67.57.Fg,67.90.+z 



The experimental realization of Bose-Einstein Conden- 
sation (BEC) in trapped alkali-atoms gases at ultralow 
temperatures offers new opportunities for studying quan- 
tum degenerate fluids j^, H, |^. The art of manipulat- 
ing these condensates, which contain thousands of atoms 
confined to microscale clouds, achieved very high level 
now. Creation of vortices in one |Q and two-component 
BEC ^ of atoms is the amazing example of this 
art. 

In this article we focus on the properties of the two- 
component BEC in the trap in which the trapping po- 
tentials for each component are displaced with respect 
to each other in the vertical direction The con- 

densate consists of the simultaneously trapped other- 
wise identical atoms of in two different hyperfine 
spin states |1 > and |2 > (|1 > and |2 > denote the 
\F = 1,771/ = -1 > and |2,1 > spin states of ^"^Rb 
atoms respectively) |@, ||, The scattering lengths of 
the states |1 > and |2 > are known to be in proportion 
Oil : ai2 : 022 — 1-03 : 1.0 : 0.97 with the average of the 
three being 55(3) A | |. 

The double condensate system was prepared from the 
single |1 > condensate by driving a two-photon transi- 
tion which transfers any desired fraction of the atoms 
to the 1 2 > state by selecting the length and amplitude 
of the two-photon pulse Q . The rotating magnetic field 
of the time-averaged orbiting potential (TOP) trap gave 
the possibility to displace the minima of the trapping po- 
tentials Vi and V2 with respect to each other. When the 
minima of the trapping potentials were not shifted, the 
|1 > atoms formed a shell about the |2 > ato ms m . This 
case has been discussed theoretically in pd|~^ If the 
minima of the trapping potentials Vi and V2 arc displaced 
from each other by a distance which is small compared to 
the size of the total condensate the resulting separation 
of the centers of mass of the condensates is much larger 



[Q . In this paper we provide for an analytical explanation 
of this result. 

In order to explore the boundary between the two con- 
densates, we begin with the analysis of their behaviour 
in the framework of the Thomas-Fermi Approximation 
(TEA), which ignores the kinetic energy terms in the 
Gross-Pitaevskii equations for the condensate wave func- 
tions . It has been shown that in the case of one com- 
ponent condensates the TEA results agree well with the 
numerical calculations for large particle numbers, except 
for a small region near the boundary of the condensate 
|l3| , In fact, even for small numbers of particles TEA 
still usually gives qualitatively correct results. 

In the dimensionless variables, the Gross-Pitaevskii 
equations for the condensates in the harmonic traps may 
be written in the form pd| : 
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l^2lVi=0; (1) 



/32V'V^ + (a:'^+/ + A2(z'-z^)2)^i- 

/ii7/.^+772/?2|V'^|Vi + ^^^^KlV^=0; (2) 



Here V'i(r) = \/ Ni/a^ip'^{r'), V'i(r) being the wave func- 
tion of the species i of a two-species condensate (i = 1,2). 
A — LOz/^- r = a_Lr', where a±_ — {h/r 



.1/2 



is the 



trapping frequency. fi[ — 2fii/hLu, where /i^ is the chemi- 
cal potential of the species i. The chemical potentials fj,i 
and yU2 are determined by the relations / = Ni. 

Ui is given by Ui = 8TTaiiNi/a±. The wave function ip[{r') 
is normalized to 1. Zq denotes the shift of the minimum 
of the trapping potential in the vertical direction. 

Equations (|l|) and (^ were obtained by minimization 
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of the energy functional of the trapped bosons given by: 
E' = \j rfV[7Vi|VV;i' + A^2/3^|VV^|2 + 
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(3) 



Here the energy of the system E is related to E' by E = 
HloE' 

In the TFA, Eqs. (|l]), (||) and (||) can be further simpli- 
fied by omitting the kinetic energy terms. In the frame- 
work of TFA the phase segregated condensates do not 
overlap, so we can neglect the last terms in Eqs. (|^), (|^) 
and (g), obtaining simple algebraic equations: 

X e(M'i-(/VA2(.' + 4f)) X 

X e(r'VA2(z'-z^)2-/4); (4) 

mirT = -U~{r'' + X'{z'~zX}) X 
X e[^,',-{r'' + X\z'-z',)')) X 
X e(r'' + A2(z'-hz^,)2-/4) . (5) 

Here Q denotes the unit step function and p'^ — x'^ -f y'^. 
If Zq — 0, from Eqs. and (j^) one can see that the 
condensate density has the ellipsoidal form. This case 
has been considered in detail in Refs. [|lO[ pi] . 

In the case of phase separation, the energy of the sys- 
tem can be written in the form |lO[ E — Ei + E2, 
where 
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M2- / dVlV'^l^ 



(6) 
(7) 



To determine the position of the boundary between 
the condensates, we use the condition of thermodynamic 
equilibrium |l5|| : the pressures exerted by both conden- 
sates must be equal: Pi — P2. Pressure is given by 
||6|: Pi = Gii\^i\'^/2, where Gu = Airh^au/mi. Using 
these equations one can obtain the equation for the phase 
boundary: 
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K — 1 



(8) 
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K = \J aii/a22, and 

j2 /^'i ~ 



(9) 



From Eqs. (||) and one can easily understand why 
the small displacement of the centers of the trapping po- 
tentials leads to the significant resulting separation of the 
centers of mass of the condensates . The basic physics 
of this amplification of the trap center difference comes 
from the two possible final configurations of the mixture 
and that the system is close to the "critical point" that 
separates the two final configurations. The two configu- 
rations are the symmetric one where one component is 
inside and the other component is outside and the asym- 
metric one||l^ [1^ where the two components are on 
opposite sides. The former configuration is favored when 
«^ = •\/ 011/022 is different from one, with the less repul- 
sive component in the middle where the density is higher. 
The asymetric configuration possess a lower interface en- 
ergy and is favored when k is close to one. We found that 
in the Thomas-Fermi approximation, when the trapping 
frequencies for the two components are the same, the 
amplification factor is proportional to 1/(k — 1). 

From Eqs. the evolution of the system upon 

increasing a may be described as follows: for a — Q con- 
densate 1 forms the shell about the ellipsoidal condensate 
2. The semiaxis of this ellipsoid is given by Eq. (||) for 
a = 0. Upon increasing a the inner ellipsoid moves up- 
wards, while external one moves down. It may be shown 
that they touch each other for the critical value of a: 

ac = \{^~ \f^\ ^ phase boundary (||) 

intersect boundaries of condensates at the points with 
coordinates: 



Xz^ 



' 1,2 
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(10) 
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which can be obtained from Eqs. (p-(||). Critical value 
ac is a function of the ratio N2/N1. 

Using normalization condition / |-i/;^(r')p cfir' = 1, one 
can determine the chemical potentials fi[ as functions of 
Ni, N2, and a. Analytical expressions for fj,'^ are different 
for a < Ur and for a > ar- In the former case one has: 
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where n° = (^^)^^^7 = = j/k. In the limit 

a — > one has the results obtained in our previous pa- 
pers for the non-displaced potential jlo[ 0. In the case 
a > ac the formulas for fi'i and fi'2 obtained after tedious 
but straightforward calculations are rather cumbersome 
and will be given elsewhere. In this article we discuss the 
results of calculations. To be specific, we will use the 
parameters corresponding to the experiments on ^^Rb 



atoms a±_ 
atoms, A 



= 2.4 X 

Vs. 



10' 



N ^ Ni + N2 ^ 0.5 X 10" 
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In Fig.l we show the density profiles of the conden- 
sates (see Eqs. (^) and (H)) as functions of the vertical 
coordinate z for r' =0 and A^i = Ni. In this case 
Q!c — 0.0047. Fig. 1(a) illustrates the experimental situ- 
ation 0: a — 0.03 (approximately 3% of the extent of 
the density distribution in the vertical direction) is larger 
than the critical value c^c^ and condensates are completely 
separated in vertical direction in accordance with the ex- 
periment [Q. In the case a < Uc the condensate is 
inside the condensate A^i . It should be noted that rather 
small shifts of the trapping potential centers with respect 
to each other produce considerable displacements of the 
condensates. The condensates in Fig.l does not overlap 
because, as was mentioned earlier, in the framework of 
TFA it is impossible to describe the overlapping of the 
condensates. 

Another interesting question is how the vortex states 
change when the minima of the trapping potentials V\ 
and V2 are displaced with respect to each other. In a 
frame rotating with the angular velocity along the z- 
axis the energy functional of the system is: 



where V';, (r) — (i")|e''^'^ is the wave function for the 
vortex excitation with angular momentum Mj. In the 
TFA, the vortex induced change in condensate density is 
negligible pO| (hydrodynamic approximation). 

In the case of the phase segregated condensate, one 
finds from Eqs. ( p^ ) and (^^ that the energy change due 
3S AE 



to presence of the vortices 
has the form l|lo|, |T1 
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and 



In the hydrodynamic limit ^[ is given by Eqs. 

(i- ■ " . 

In the case a = critical velocities as functions of ratio 
N2/N1 have been calculated in |l^. It was shown that 
for all values of N2/N1, the critical velocity JIat^ of the 
inner condensate is lower than the critical velocity SljVi 
of the outer one. So upon increasing fi, a vortex will 
appear first in the external condensate. However, if for 
given one shifts the centers of the trapping potentials 



with respect to each other in the vertical direction, the 
inner condensate floats to the surface. In this case one 
can expect that the critical velocities of the condensates 
become closer and even can be equal for some values of 
a and N2/N1. 

As in the case of the chemical potentials, the expres- 
sions for the critical velocities have different analytical 
forms for a < ac and for a > ac, the latter being rather 
cambersome. For a < ac the critical velocities are given 
by: 



1(3/2) 
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(16) 



(17) 



Again, in the limit a — *■ one has the results obtained in 
our previous paper for the non-displaced potential | |To[ |. 

Figure 2 shows the behaviour of critical velocities as 
functions of a for different values of N2/N1. Dashed lines 
correspond to the inner condensate, solid lines - to outer 
one. From Fig. 2(c) one can see that the critical velocities 
really can intersect. Physically this means that there is 
a possibility to exchange the vortex states between con- 
densates by shifting the centers of the trapping potentials 
with respect to each other for fixed angular velocities. 

In summary, we investigated the behaviour of simulta- 
neously trapped codensates consisting of atoms in 
two different hyperfine states. It is shown that the small 
shift of the minima of the trapping potentials with re- 
spect to each other leads to the considerable displacement 
of the centers of mass of the condensates in agreement 
with the experiment . It is also shown that the critical 
angular velocities of the vortex states of the condensates 
drastically depend on the shift and relative number of 
particles in the condensates. The predicted exchange of 
the vortex states between the condensates as a function 
of the shift remains to be studied experimentally. 
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Figures captions 

Fig. 1. Density profiles of the condensates as functions 
of the vertical coordinate z for A^i — N2. Figure (a) 
corresponds to a > ac and figures (b)-(c) - to a < ac- 
Sohd lines correspond to the 1 1 > atoms and dashed lines 
- to the 1 2 > atoms. 

Fig. 2. Critical velocities of outer condensate VIn^/u) 
and the inner condensate Qn^/u as functions of a for 
different values of N2/N1. Dashed lines correspond to 
Q.f4^/uj and soHd hues - to OjVi/t^- 




O lA O lA O lA OO O lA O lA O lA OO O lA O lA O lA OO 
CO CM CM 1- 1- O O CO CM CM 1- 1- O O CO CM CM -i- -i- O O 

OOOOOOO OOOOOOO OOOOOOO 



